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Abstract 

In 1958 Jeffreys proposed a power law of creep, generalizing the 
logarithmic law earlier introduced by Lomnitz, to broaden the 
geophysical applications to fluid-like materials including igneous rocks. 
This generalized law, however, can be applied also to solid-like 
viscoelastic materials. We revisit the Jeffreys-Lomnitz law of creep 
by allowing its power law exponent a, usually limited to the range 
< a < 1 to all negative values. This is consistent with the 
linear theory of viscoelasticity because the creep function still remains 
a Bernstein function, that is positive with a completely monotone 
derivative, with a related spectrum of retardation times. The complete 
range a < 1 yields a continuous transition from a Hooke elastic solid 
with no creep (a— )■— oo) to a Maxwell fluid with linear creep (a = l) 
passing through the Lomnitz viscoelastic body with logarithmic 
creep (a= 0), which separates solid-like from fluid-like behaviors. 
Furthermore, we numerically compute the relaxation modulus and 
provide the analytical expression of the spectrum of retardation times 
corresponding to the Jeffreys-Lomnitz creep law extended to all a < 1 . 
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Transform. 
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1 Introduction 



According to the linear uni-axial theory of viscoelasticity, a material body can 
be viewed as a linear system where either stress a{t) or strain e{t) have the 
role of excitation function (input) and response function (output). From both 
the experimental and theoretical viewpoints, a central role is played by the 
creep and relaxation tests, where the inputs are a step-wise stress and strain, 
respectively. Introducing the Heaviside step function 6(t), in a creep test 
a{t) = cro6(t), conversely in a relaxation test e{t) = eoO(t), with constants 
(Jo and 69. The corresponding outputs are characterized by the so-called time- 
dependent material functions, the creep compliance J{t) = e{t)/aQ and the 
relaxation modulus G{t) = (T(t)/eo, respectively. From experimental evidence, 
J{t) is non- decreasing and non-negative, while G{t) is non-increasing and 
non-negative. 

In Earth rheology, transient creep is often described by the Jeffreys- 
Lomnitz law. This law of creep is defined in terms of four parameters as 
follows 

J(t) = Jo[l + #(t)], t>0, (1) 

where Jo is the unrelaxed compliance, g is a positive dimensionless material 
constant, and 

(i.l)"-i 

^(t) = A ^ 0<a<l, t>0. (2) 

a 

Here, Tq is a characteristic time and the exponent a essentially determines 
the rheological behaviour of the creep law. By its own definition, ilj{t) is 
a positive increasing function of time that we refer to as the dimensionless 
Jeffreys- Lomnitz creep law. 

Creep law ([1]) was proposed by Sir Harold Jeffreys in 1958 [21] (see also [25] 
and [261 122]), to generalize, via the parameter a, the logarithmic law 

^(t) = log (^1 + -^^ , t>0. (3) 

introduced by Lomnitz in 1956 [29] to describe flow in igneous rocks. This 
law was also employed by Lomnitz to account for the damping of the Earth's 
free nutation (Chandler wobble) and of seismic S-waves, see Lomnitz (1957, 

1962) [sniEi]. 
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As a matter of fact, for a = 0, the Jeffreys law reduces (in the hmit) to the 
Lomnitz logarithmic law; while for a = 1, it reduces to the simple linear 
law of the Maxwell viscoelastic body. According to Jeffreys, Eq. ([2]) better 
fits the data of creep and dissipation in the Earth for seismological purposes. 
The Jeffreys-Lomnitz law is indeed well recognized in rheology of rocks, in 
view of its property of interpolating creep data between a logarithmic and a 
linear law, see e.g. Jeffreys and Ranalli 

Although in the original papers the exponent a was not explicitly limited 
to positive values, the applications with a < were only consequently 
introduced in Earth rheology. We note that the late Professor Ellis Strick 
in his 1984 paper |33] had already extended the Jeffreys-Lomnitz law of 
creep allowing the parameter a in the range — !<«<+!, by introducing 
parameter s:=l — a;(0<s<2). Furthermore, he was also interested in 
the representation of the extended Jeffreys-Lomnitz law of creep in terms 
of a suitable ladder network of springs and dashpots. In his work, Strick 
was motivated by some experimental observations suggesting negative values 
of the exponent a. At the time we started editing the present paper, we 
were not aware of any reaction of the geophysical community to his results. 
However, later we noted some published works, see e.g. Crough & Burford 
(1977) [S], Spencer (1981) [H], Wesson (1988) Darby & Smith (1990) 
[H], where the Jeffreys-Lomnitz law of creep was applied with a < to fit 
experimental geophysical data. In these papers, however, we did not find any 
attempt to characterize the viscoelastic properties of the Jeffreys-Lomnitz 
law in the whole range — cxd < a < 1. In particular, the relaxation modulus 
and the retardation spectrum associated to this creep law have never been 
investigated up to now. Of course, in the trivial case a = 1, we recover the 
Maxwell body for which the retardation spectrum vanishes and the relaxation 
modulus is exponentially decreasing. 

The purpose of this paper is to characterize the Jeffreys-Lomnitz law of 
creep consistently with the linear theory of viscoelasticity. In order to 
meet this goal, we find it convenient for the readers to recall in Section 
2 the essential mathematical notions of linear viscoelasticity necessary to 
understand our analysis. Then, in Section 3, we consider this law of creep 
allowing the parameter a in the whole range — oo < a < 1. Here, for 
some selected values of a in the range —2 < a < 1, we show plots of the 
analytically determined creep function and of the corresponding relaxation 
function, numerically computed by solving an integral Volterra equation 



3 



of the second kind. In Section 4, we analytically derive the spectrum of 
retardation times corresponding to the extended creep law by using the tools 
of Laplace transforms, and we exhibit plots for the above selected values of a. 
In Section 5, we draw our conclusions and final remarks. We have assumed 
that the reader is familiar with the main properties and tables of Laplace 
transforms so we do not recall them. However, we devote the Appendix to 
the Post-Widder formula of Laplace inversion on the real axis, which yields 
a complementary approach to our derivation of the retardation spectrum. 

2 Essentials of linear viscoelasticity 

The basic principles of linear viscoelasticity are illustrated in well known 
treatises [SHI IQ [HI [3 |36l gGl |28]. Here we mostly follow the notation of 
Chapter 2 of the recent book of Mainardi [33], where the approaches by 
Gross [20], Pipkin 06] and Tschoegl [16] have been mainly adopted and 
revisited. Without loss of generality, the one-dimensional theory will be 
employed throughout. 

Assuming causal and differentiable histories of stress and strain (i.e., aif) 
and e(t) vanish for all times t < and are differentiable for t > 0), the 
time-dependent material functions formerly introduced, J if) and G{t\ can 
be used to express the general stress-strain relationship as 

tit) = (t(0+) J{t) + f Jit- t) (T(r) dr = (t(0+) J(t) + J{t) * a{t), (4) 
Jo 

or, equivalently 

a{t) = e(0+) G{t) + I G{t- r) e(r) dr = e(0+) G{t) + G{t) * e{t), (5) 
Jo 

where dots and * denote time-differentiation and time-convolution, respec- 
tively. 

The limiting values of the time- dependent material functions for t — )■ O"*" 
and t — )■ +00 describe the instantaneous and equilibrium behaviour of the 
viscoelastic body, respectively. The instantaneous Jo and the equilibrium 
comphance Joo are defined by Jo = -'^(O^) and Joo = J(+oo), while the 
instantaneous and equilibrium modulus are Go = G'(0"^) and Goo = G(+oo), 
respectively. 
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Taking the Laplace transform of Eqs. and the following reciprocity 
relation is obtained 



sJ(s) = — i— or Jis)G{s) = \, (6) 



sG{s) 



s 



where s is the Laplace complex variable and the tilde denotes transformed 
quantified. Applying the limiting theorems of the Laplace transform in the 
first of Eqs. provides 

Jo = and = , (7) 

which allow us for a classification of the viscoelastic bodies into four types, 
according to the values of moduli and compliances (Tabled]). 



Table 1: Classification of the four types of viscoelastic bodies according to 
values of Jo, Joo,Go and Goo (conventionally, 1/0 = oo). 



Type 


Jo 


J DO 


Go 


Goo 


I 


> 


< oo 


< oo 


> 


II 


> 


= oo 


< oo 


= 


III 


= 


< oo 


= oo 


> 


IV 


= 


= oo 


= oo 


= 



Type I bodies exhibit both instantaneous and equilibrium elasticity, so they 
behave similarly to an elastic (Hooke) body for sufficiently short and long 
times. Bodies of type II and IV show a complete stress relaxation (at constant 
strain) since Goo = and an infinite strain creep (at constant stress) since 
Joo = oo,, so they do not exhibit equilibrium elasticity. Finally, bodies of type 
III and IV do not present instantaneous elasticity since Jq = (Gq = oo), 
hence they will not be considered in this study. By the Laplace inversion 

^For the most relevant properties of the Laplace transform, we refer, e.g,. to the 
Appendix A3 (Eq. 32) of the well-known treatise by Tschoegl, see pp. 560-570. In 
the following we use the notation f{t) ~ /(s) to denote the juxtaposition of an original 
function f{t) with its Laplace transform /(s). 
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of the second of Eqs. (jS]), we obtain the interrelation between the material 
functions in the time domain 



J{t) * Git) := [ J{t-t')G{t')dt' 



= t, t>0, 



(8) 



Jo 



which will be useful in the following to obtain, by numerical methods, G{t) 
from the knowledge of J{t) in the case of the (generalized) Jeffreys-Lomnitz 
law of creep. In fact, if Jq > (types I and II), Eq. ([8]) can be re-written as 
a Volterra integral equation of the second kind 



with G{t) unknown. Of course, a similar equation can be established for 
J{t), provided that G(t) is known and Gq < oo (types I and II). 
Another relevant consequence of Eq. (|H]) is 



where the equality holds in the limiting cases t — )■ 0^ and t — ?► +oo. The 
results in Eqs. (E]), ( fTOl) are discussed in [36] and recently revisited in [33] . 

We remark the restrictive conditions related to the monotonicity of the 
time-dependent material functions. According to Gross pO], their general 
expressions consistent with experimental observations are as follows: 



where the non-negative functions -Re(T) and Raij) denote the retardation 
spectrum and the relaxation spectrum^ respectively. For the classical 
mechanical models, obtained by series and parallel combinations of a finite 
number of springs and dashpots [33], the spectra are discrete. 

It is convenient to consider separately, for t > 0, those terms of (fTT]) deriving 
from a continuous spectrum 





Jit) Git) < 1 , t > , 



(10) 




(11) 




(12) 
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where Jr{t) (the creep function with spectrum) is a non- decreasing, non- 
negative function with hmiting values Jt{0~^) = and J^(+oo) < oo, 
whereas Gr{t) (the relaxation function with spectrum) is a non-increasing, 
non-negative function, with ^^(0"*") < oo and Gr{+oo) = 0. The integrabihty 
on of Rei^r) and Rcf{t) is ensured if Jr(+oo) < oo. 

We also note that, consistently with the spectral representations (fT2|) . 
we have (— l)" Jr"'^(t) < and {—1)'"-G^r\t) > 0, where the superscript 
denotes the n-th time derivative. From a mathematical standpoint, these 
conditions are equivalent to require that Jrit) and Gr{t) are Bernstein and 
complete monotone functions, respectively. For details, we refer to specialized 
mathematical treatises. Berg and Forst [S], Gripenpeg et al. [12], and 
Schilling et al. [IQ]. These properties have been investigated by several 
authors, including Molinari [33] , Del Piero and Deseri [10] and, more recently, 
by Hanyga [21]. 

The determination of the time-spectral functions from the knowledge of the 
time-dependent material functions is a problem that can be formally solved 
through a method, outlined by Gross [20] , based on Laplace transform pairs. 
For this purpose we introduce the frequency- spectral functions as 

S.W:^^. (13) 

where 7 = 1/r denotes a retardation or relaxation frequency. Noting that 
Reij) dr = Sei'y) d'j and R^ij) dr = Sa{'y)dj, time differentiation of ( 1T2|1 
yields 



Jo r Jq 



-tl 



d-f, 



Gr (t) = I ^ e-^/^ dr= (7) e-*^^ d^ 



(14) 



showing that 75*^(7) and 7 6*0- (7) can be viewed as the inverse Laplace 
transforms of Jrif) and —Grit) , respectively, where t is now considered 
the Laplace transform variable instead of the usual s. Thus, adopting the 
connective symbol -r- for Laplace transform pairs where in the LHS we put 



7 



the original function and in the R.H.S. its Laplace transform, we have 



+7^e(7) 



7 



-iSaKl) = ^ Grit) . 

7 



(15) 



Consequently, when Jr(t) and Gr(t) are known analytically, the correspond- 
ing frequency-spectral functions can be derived by standard methods for 
the inversion of Laplace transforms; then, from Eq. f ll3p . the time-spectral 
functions are easily obtained. 



3 The extended JefFreys-Lomnitz laws of 
creep and relaxation 

In this Section we revisit the Jeffreys-Lomnitz law of creep by extending the 
range of variability of the exponent a to any negative value in Eq. ([2]). In 
the expression for the dimensionless Jeffreys-Lomnitz creep law, it is now 
convenient to separately consider four cases: 



' t/To, 

(1 + t/ror-l 



a 



1 



a 



t>o, i){t) = < 



log(l + t/ro) , 
l-(l + t/ro) 



< a < 1 



a 



a <0. 



(16) 



\a\ 



The behaviour of ip{t) as a function of the dimensionless time t/ro is 
illustrated in Figure [H for some values of a in the range —2 < a < 1, 
adopting a logarithmic and a linear time axis in the top and in the bottom 
frames, respectively. From Eq. f ll6p . the different asymptotic behaviour of 
ip{t) as t — 7- oo in the cases < a < 1 and a < is apparent. In the former 
case, ilj{t) diverges whereas in the latter case it increases up to the finite 
value l/\a\. 
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Figure 1: The creep function ipi^t) for a = 1, 0.5, 0, —0.5, —1, —1.5, —2 versus 
dimensionless time t/ro, adopting a linear scale (top) and a logarithmic scale 
(bottom) . 
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In consideration of the asymptotic behaviour of ip{t), from the classification 
of the types of viscoelasticity of Table [H we recognize that the law provided 
by Eq. f|T6|) . along with ([1]), describes the behavior of a viscoelastic body that 
shows instantaneous elasticity. For a < and < a < 1, it falls in types I 
and II, respectively. 

Since the extended Jeffreys-Lomnitz law accounts for instantaneous elasticity, 
it is possible to write the relaxation modulus as follows 

G{t) = Go<P{t) , t>0, (17) 

where 0<Go = l/^o<oo and (pit) is a positive decreasing function of time 
that we refer to as the dimensionless relaxation function. 

Using Eqs. ([1]) and (fT7|) in (Q, the following interrelation in the Laplace 
domain is found between the dimensionless creep and relaxation functions: 

s0(s) = or 0(s) = (18) 

1 + qs ip^s) 1 + qsipis) 

where 0(s) and ipi^s) are the Laplace transforms of (j)(t) and '?/'(t), respectively. 
From the limiting theorems of Laplace transforms applied to the first of Eqs. 
(fT8|) . we easily obtain the following: 

f 0, < a < 1, 

= ^(+00) = I (19) 

[ l + 9/|ar 

However, the analytical derivation of (f){t) from inverting the Laplace 
transforms in Eqs. f|T8l) appears as a difficult (if not prohibitive) task even 
if ip{s) can be expressed in terms of a transcendental function related to 
incomplete Gamma function, see [Hj. In order to obtain the relaxation 
function, we numerically solve the Volterra integral equation of the second 
kind ([9]) with J{t) appropriate for the extended Jeffreys-Lomnitz law. After 
straightforward algebra, the Volterra equation reads: 

0(t) = l-g^*(^l + ^^ <Pit')dt'. (20) 

From now on, for the sake of simplicity, we will assume q = 1- The numerical 
method is implemented in the Fortran routine voltra.for, available from 
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Numerical Recipes [2Z]. We note that an alternative semi-analytical method 
oriented to viscoelasticity has been first introduced by Hopkins and Hamming 
in the late 1950s |22l |23], and recently revisited by Lin in his book [28]. We 
also note that, in the limit a — >■ — cxo, = 1 and the elastic Hooke model 
is recovered. 

In Figure we show the relaxation function versus the dimensionless 
time t/ro, for some values of a in the interval —2 < « < 1, adopting a 
logarithmic and a linear time axis in the top and in the bottom frames, 
respectively. 

We conclude this section with a discussion of the asymptotic representations 
of and as t/ro — )■ 0+ and t/ro — )■ +00, which provide simple 
analytical estimates of the corresponding functions for sufficiently small and 
large times, respectively. These asymptotic results, indeed trivial for V'l^), are 
more relevant for for which we only dispose of a numerical evaluation. 
From the asymptotic analysis we exclude the trivial case a = 1, for which 
we easily recover the creep and relaxation laws of the Maxwell body, valid 
for any time: 

^{t) = t/ro , = exp(-t/ro) t > 0. (21) 

For the function ipit), the case t/ro ~^ 0+ can be easily studied recalling the 
binomial series representation 

which is convergent for < t/ro < 1. Hence, we obtain the straightforward 
asymptotic result 

t/ro ^ 0+ ^ tpit) ~ t/ro , « < 1 . (23) 

In the case I/tq +00, starting from the definition given by Eq. flTB]) . we 
have: 



i/TQ +00 ^> il){t) ~ < 



' (t/ro)"/a, < a < 1, 

log(t/ro), a = 0, (24) 

^ (1 - (t/ro)-H) /|«| , «<0. 
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0.0 

1 2 3 4 5 

t/To 




LOGio(t/To) 



Figure 2: The relaxation function 4>{t) for a = 1, 0.5, 0, —0.5, —1, —1.5, —2 
versus dimensionless time t/ro, adopting a linear scale (top) and a logarithmic 
scale (bottom). 
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For the function (p(t), the asymptotic results for t/rg — )■ 0+ and t/vQ — )■ +00 
can be obtained by applying the Tauberian theorems in the R.H.S. part of 
Eq. f lTSj) (with g = 1) of Laplace transform. Because we prefer to avoid the 
exact expression of ipi^) terms of transcendental functions, we will use the 
Tauberian theorems in a simplified form. This means to invert for s ^ 00 
the Laplace transform of ?/'(t) as t — )■ 0, and conversely for s — ?■ the Laplace 
transform of '?/'(t) as t — t- 00. 

As far the behaviour as t/ro — )■ 0^ is concerned, from Eqs. (123|) and (fT8|) we 
have: 



s — )■ 00 



is) 



1/s 



1/s 



1/s -1/s' 



so that 



t/ro ^ 0^ 



~ 1 - t/ro , a <1. 



(25) 



As far the behaviour as t/ro ~^ +00 is concerned, from Eqs. ([21]) and f llSp 
we obtain: 



f r(a)/s 



s^O 



a+l 



i/j{s) ~ < 



-{C + \ogs)/s . 



< a < 1 



a = 



^ l/(|a|s) + r(-|a|)/s-l"l+S a<0 



and 



1/s 



[s] ~ < 



l + r(a)/s° r(a)si- 
1/s 1 



1 — C — logs slog(l/s) 
^ l + l/\a\+T{-\a\)/s-\"\ ' 



< a < 1 



a = 



a < , 



where for the case a = the constant C = — r'(l) = 0.577215... denotes the 
so called Euler-Mascheroni constant. 
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Henceforth 



sin(7r«) 



t/ro -> +00 ^> (f)(t) ~ < 



a 



< a < 1, 



a = 0, 



log(t/ro) ' 

T-j^fl + T^(Vro)-l"lY «<0 
{ 1 + l/\a\ \ 1 + |a| J 



(26) 



where the Laplace inversion of l/(s log(l/s)) was obtained via the Karamata- 
Tauberian theory of slow varying functions in the paper by Mainardi et 
al. [32]. For details on slow varying functions and on Karamata-Tauberian 
theorems, see e.g., Feller (1971) [T7] . 



4 Retardation spectrum for the extended 
JefFreys-Lomnitz law of creep 

Following the method outlined in Section 2, see the top Eqs. in ( IT4|) or 
(IT^ . we must consider the derivative of the dimensionless creep function in 
the extended Jefferys-Lomnitz law as the Laplace transform (with Laplace 
parameter ^ = t/rg) of an "original function" in the variable 7 = To/t. 
For sake of simplicity, from now on, we will scale times and frequencies by 
assuming tq = 1, and also put the material constants Jq and q equal to 1. 

In order to obtain the required time retardation spectrum we first derive 
the corresponding frequency spectrum for a < 1 by proving the Laplace 
transform pair: 

^S.(7) = f^^=.K + ir- = | = ^.tt). a<l. (27) 

The prove is based on the well-known integral representation of the Gamma 
function, 

poo 

T{z):= e-"" u^'-^Uu , Re{4 > . 
Jo 
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In fact 



1 r°° p-7 ((: _|_ 1 foo 



r(l - a) Jo 7" r(l-a 

where we have used the change of variable m = + 1)7 and the integral 
representation of r(l — a). Then, from the first expression in Eq. (!T3|l we 
have 

/?,(l/7)=7'5.(7) = :p7T^^, «<1, 

i (1 — a) 7" ^ 

so that, noting r = I/7, we derive the required analytical form for the 
retardation spectrum. 

For the limiting values of a = 1 (Maxwell body) and a — t- —00 (Hooke body) 
we recover Rei^) = due to the vanishing of the factor l/r(l — a). 

Accounting the previous notes on the asymptotic behavior of iplt) we now 
explore the integrability of such spectrum on the interval < r < 00. Using 
directly Eq. ( I28ll . we have 



Tr^W 1 r^"" . / ifO<a<l, 

io^^^"^ rxI^yo^ '^lVH if-oo<«<0. (29) 

This is obtained by changing the variable r into u = 1/r, recalling the integral 
representation of the Gamma function and noting r{—a) = — r(l — a)/a. 

Of course the above result is consistent with the asymptotic behaviour of our 
creep law (fT6l) as t/ro +00. 

In Figure 131 we show the retardation spectrum -Re(r) versus the dimensionless 
retardation time t/tq for some values of a G [— 2, +1], adopting in abscissa 
a logarithmic scale. 

We note that, for a = (Lomnitz body), the maximum of the spectrum 
occurs just for t/tq = 1. Furthermore, for a < 0, as r — the maximum 
increases with |a| up to infinity in such a way that identically vanishes in 
the open interval r > 0. 
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0.75 




LOGioC-c/To) 

Figure 3: The retardation spectrum Re{T) versus dimensionless time t/tq, 
for a = 1, 0.5, 0, —0.5, —1, —1.5, —2, adopting a logarithmic scale. 



5 Concluding remarks 

By revisiting the Jeffreys-Lomnitz law of creep, we have described the 
behaviour of a general viscoelastic body that, according to the value of a 
single parameter a ranging from 1 to — oo, shows a transition from a Maxwell 
to a Hooke body. This continuous transition is well shown in Figured! where 
the creep law ilj{t) is displayed for some values of a in the range —2 < a < 1 
for more rheo logical interest. We have numerically computed the relaxation 
modulus (f){t) corresponding to our creep law by solving a Volterra integral 
equation of the second kind, see Figure O Furthermore, we have provided the 
analytical expression of the spectrum of retardation times, which is shown in 
Figure [31 

We are not aware of any previous attempt in the literature to compute the 
retardation spectrum and the relaxation law for the Jeffreys-Lomnitz law 
of creep, so we think to have filled a gap in this field. However, we have 
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not numerically evaluated the relaxation spectrum leaving this problem to 
interested readers who can chose the most suitable method to approximate 
the relaxation spectrum from the numerical data of the relaxation modulus. 
On this respect there exists a huge literature, see e.g. [U |2l |3l HJ [TTl [121 [131 

We also note that in the recent book by Schilling et al. [IQ], where 
a rigorous mathematical description of Bernstein functions is found, the 
authors have considered a function similar to our extended law of creep in 
two examples (see No. 2,3 at pp. 218-219, Chapter 15) without the factor 
l/a, corresponding to our parameter < a < 1 and — 1 < a < 0. For these 
cases they have provided the so-called Levy density that, as a matter of fact, 
coincides with our frequency spectral function S^i'j) and therefore with our 
retardation spectrum Re{T) by putting 7 = 1/r. Therefore, in this work, we 
have extended their results to the cases a = and a < — 1 by adopting our 
notation. 

Finally, we hope that the results obtained for the creep and relaxation laws 
investigated in this paper may be useful for fitting experimental data of creep 
and/or relaxation responses in rheology of real materials. 
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Appendix: The Post-Widder formula for the 
retardation spectrum 

The Post-Widder formula provides the original function f{t) from its Laplace 
transform /(s), known with all its derivatives on the real semi-axis of the 



17 



complex Laplace plane, via a limit of infinite sequence as 

.l)n 



fit) = lim 



s=n/t 



(Al) 



where is the n-th derivative of / with respect to the Laplace variable 
s. This is the case when the Laplace transform is proved to be an analytic 
function on the right-half s-plane. For other details we refer e.g. to |46j. 

The purpose of this Appendix is to prove the validity of Eq. (!27|) by using 
(in a suitable way) the Post-Widder formula. In our case, denoting in Eq. 
(A.l) t by 7 and s by C,, we must verify that 
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with /(O := + 1)""^ and a < I. Now 
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so we get 
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We easily recognize from Stirling asymptotic formula that for n — )■ oo 
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we finally get 
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in view of the Neper limit 
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